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^^ Abstract This chapter reviews the occurrence of quantised vortices in polariton flu- 

I— I ids, primarily when polaritons are driven in the optical parametric oscillator (OPO) 

regime. We first review the OPO physics, together with both its analytical and nu- 
merical modelling, the latter being necessary for the description of finite size sys- 
tems. Pattern formation is typical in systems driven away from equilibrium. Simi- 
larly, we find that uniform OPO solutions can be unstable to the spontaneous for- 
mation of quantised vortices. However, metastable vortices can only be injected ex- 
ternally into an otherwise stable symmetric state, and their persistence is due to the 
OPO superfluid properties. We discuss how the currents charactering an OPO play 
a crucial role in the occurrence and dynamics of both metastable and spontaneous 
vortices. 



o 

O 1 Introduction 



Quantised vortices are topological defects occurring in macroscopically coherent 
systems, and as such have been broadly studied in several area of physics. Their 
QQ existence was first predicted in superfluids |IT||2), and later in coherent waves |3J. 

■^ Nowadays, quantised vortices have been the subject of extensive research across 

■^ several areas of physics and have been observed in type-II superconductors, '^He, 

r — ultracold atomic gases, non-linear optical media (for a review see, e.g., [4]) and 

t^ very recently microcavity polaritons ||5] |6l |7] |8l |9l [TO] [TT] [T2|, the coherent strong 

^~~^ mixing of a quantum well exciton with a cavity photon. 
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This chapter reviews the occurrence of quantised vortices in polariton fluids, pri- 
marily when polaritons are driven in the optical parametric oscillator (OPO) regime. 
The interest in this area of research is manifold. To start with, the search for conden- 
sation in solid state excitonic systems has been arduous and lasted more than two 
decades: Unambiguous evidence for condensation has been reported for microcav- 
ity polaritons for the first time in 2006 1 1 3 1 . These results have been followed by a 
wealth of experimental and theoretical advances on aspects related to macroscopic 
coherence, condensation, superfluidity, quantum hydrodynamics, pattern formation, 
just to mention few (for a review, see Ref. |14|). Two different schemes of inject- 
ing polaritons and spontaneously generating a macroscopically coherent state can 
be employed: (i) non-resonant pumping, and (ii) parametric drive in the optical- 
parametric -oscillator (OPO) regime. What both condensates have in common is 
the phenomenon of spontaneous phase symmetry breaking (and the consequent 
appearance of a Goldstone mode), and the non-equilibrium ingredient. However, 
the way polaritons are pumped has strong effects on the type of condensed regime 
that can be reached. In both regimes (i) and (ii), the quest for superfluid behaviour 
has been and is being widely investigated. How it has been recently discussed in 
Refs. lfT5][T6l fl7|, one of the aspects that makes condensed polariton systems novel 
compared to known superfluids at thermal equilibrium, is that now all the paradig- 
matic definitions of a superfluid, such as the appearance of quantised vortices, the 
Landau criterion, the existence of metastable persistent flow, the occurrence of soli- 
tary waves, have to be singularly examined and might in general be fundamentally 
different from the equilibrium case. Several of these popular topics are examined 
in other chapters of this book cross-refer to : A. Bramati and A. 
Amo, B. Deveaud, D. Krizanovskii and M. Skolnick, F. Laussy, 
G. Malpuech, M. Wouters and V. Savona, D. Snoke, Y. Yamamoto. 

Resonantly pumped polaritons in the OPO regime ITSl IT9l have been recently 
shown to exhibit a new form of non-equilibrium superfluidity ll20l l8l. Polaritons 
continuously injected into the pump state, undergo coherent stimulated scattering 
into the signal and idler states. Superfluidity has been tested through as frictionless 
flow II20I of a travelling signal triggered by an additional pulsed probe laser (the 
TOPO regime). In addition, the study of quantised vortices imprinted using pulsed 
Laguerre-Gauss laser fields has attracted noticeable interest both experimentally |8 1 
and theoretically 112111221 |23ll24l . providing a diagnostic for superfluid properties of 
such a non-equilibrium system. In particular, vorticity has been shown to persist not 
only in absence of the rotating drive, but also longer than the gain induced by the 
probe, and therefore to be transferred to the OPO signal, demonstrating metastability 
of quantised vortices and persistence of currents Il8l l22l . 

The chapter is arranged as follows: after a very short introduction to microcavity 
polaritons in Sec.l2] we describe the optical parametric oscillator regime in Sec. [3] 
stressing the analogies and differences with an equilibrium weakly interacting Bose- 



Einstein condensate (Sec. 3.2.1 1 and the numerical modelling that is necessary to 



use for finite size pumps (Sec. 3.3 1. In Sec. 3.5 the occurrence of spontaneous stable 



vortices in OPO is described for clean cavities, while the case of disordered cavities 



is studied at the end of Sec. 3.3 Next we describe in general terms the role of 
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adding a pulsed Gaussian probe to the OPO regime (the so called TOPO regime) in 
Sec. [4] while metastable vortices triggered by a Laguerre-Gauss probe are discussed 
in Sec. l5] Here, in Sec. 5.2 we also describe the onset and dynamics of vortex- 



antivortex pairs. Stability of multiply quantised vortices is analysed in Sec. l6] and 
finally we mention the occurrence of vortices in polariton fluids in other regimes 
than OPO in Sec. El 



2 A very short introduction to microcavity polaritons 

Before focusing on the main topic of this review, we give here a very short in- 
troduction to microcavity polaritons in order to fix the notation for later on. A more 
complete introduction can be found in several review articles ll25ll26ll27ll28ll29l[T4l 
and books ll30l ISTl l32l [33 1 on microcavity polaritons. 

Microcavity polaritons are the normal modes resulting from the strong coupling 
between quantum well (QW) excitons and cavity photons. In semiconductor micro- 
cavities, the mirrors employed to confine light are distributed Bragg reflectors, i.e., 
alternating quarter wavelength thick layers of dielectrics with different refractive 
indices. Between the Bragg reflectors, the cavity light forms a standing wave pat- 
tern of confined radiation, which can be described by an approximatively quadratic 
dispersion, (Oc{k) = (0^ + k^/{2mc) (from here onwards we fix /i = 1). Excitons 
are the hydrogenic bound states of a conductance band electron and a valence band 
hole, therefore their mass is much larger than the cavity photon mass (typically 
mc ~ lO^^iUe, where ifig is the free electron mass). For this reason, exciton dis- 
persion can be neglected, (Ox{k) = (0^. In microcavities, one or multiple QWs are 
grown in between the mirrors, so that excitons are at the antinodes of the confined 
light, giving rise to strong coupling. In addition, cavity mirrors are built with a 
wedge, so as to change the detuning between the normal incidence energy of the 
cavity field and the exciton one, 5 = 0)^ — O^. Typical parameter values for a GaAs- 
based microcavity are listed in Tab. [T] 

The polariton normal modes can be found by solving the coupled Schrodinger 
equations for exciton and photon fields, y/x.c — Vx.c(r,f) 

id, (":') = Ho M i/o = f 1-^ „ . ^^(' ..„ 1 , (1) 



WcJ \WcJ V ^«/2 «c(-'V)-/K-c 

where Qr is the Rabi splitting and Kx.c are the decay rates of exciton and photon. 
For an ideal cavity, Kx.c = 0, the eigenstates of this equations in momentum space, 
V/A:,c(r,f) = e"^''Lk^'^'^Wx.c.ky ai"e the lower (LP) and upper polaritons (UP) 
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Fig. 1 Lower (LP) and upper polariton (UP) energy dispersions (solid black), together with the 
dispersions of the photon (C) and exciton (X) fields (red dashed) as a function of either the wave- 



vector k [flm ] or the emission angle (p [degree] for m^ 
1.5280 eV, and a detuning 5 = 1 meV. 
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with an energy dispersion given by (see Fig.fTll: 
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At zero detuning (5 = 0) and normal incidence {k = 0) polaritons are exactly half- 
light and half-matter quasi-particles (cos^ 0o — 0.5 — sin^ do). The value of the mo- 
mentum k of polaritons inside the cavity is related to the emission angle (p (with 
respect to normal incidence) of photons outside the cavity by ck = a)Lp(k)sin^. 
Thanks to this property, microcavity polaritons can be directly excited by a laser 
field and detected via reflection, transmission or photoluminescence measurements. 
In Fig. [T] the energy dispersion of the lower and upper polariton are plotted as a 
function of both wave-vector, k [/xm^'], or the emission angle, ^ [degree], for typi- 
cal values of microcavity parameters. 
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2.1 Exciton-exciton and exciton-photon interaction 



A fundamental property of polaritons is their non-linear behaviour inherited from 
the exciton-exciton interaction and the saturation of the exciton-photon coupling. In 
this review, we treat excitons as bosonic particles, therefore the effective exciton- 
exciton interaction can be written as 






XX 



1 

24 



E ^?V4,k+qV4,k'-q'fe,kV0f,k' 



k,k',q 



where the effective interaction potential Vc/ can be determined starting from the mi- 
croscopic electron-hole Hamiltonian |34, 35 1. The typical wave-vectors involved 
in the physics described by this review are much smaller than the inverse exci- 
ton Bohr radius, q ^ a^', where ax = e/(2/ze^) is the two-dimensional exci- 
ton Bohr radius, e the dielectric constant, and /i^' = m^' +"V '^he electron- 
hole reduced mass. In this limit, it can be shown ll34l that the momentum depen- 
dence of Vq can be neglected, thus approximating it with a contact interaction, 
Vg ^ gx = 6e^ax/e = 6^y;f a|, where ^y^ = e^/{eax) = l/{2pLa\) is the ex- 
citon Rydberg. Typically, for GaAs quantum wells (see Tab. [Til, £ = 13, ox — 7 nm, 
and S^yx — 17 meV , therefore gx — 0.005 meV(/im)^. We will see, however, that 
the exact value of the coupling constant gx has no relevance for the mean-field dy- 
namics we are going to describe, i.e., gx can be rescaled to 1. 



QW 


cavity 


£-13 


(B^~co^~ 1.53 eV 


Me = 0.063m° 


de [-10,10]meV 


Mh = 0.3m? 


mc = 2.3 X IQ-^nP, 


ax^VOi 


Ic = 0.868 iim 


Myx ^ 17 meV 


Qr ~ 4.4 meV 


Kx ~ /xeV 


Kc = 0.1 meV 



Table 1 Characteristic parameters of a GaAs-based microcavity, divided between the parameters 
of the quantum well (left) and those describing the microcavity (right). Here, (c = \/l/(mc^fi) 
is a characteristic length for the cavity photons. The photon decay rate Kc refers to a cavity mirror 
with typically 25 bottom pairs and 15 lower pairs (see, e.g., Ref. 1 8 1). 



The composite nature of excitons, as a bound state of an electron and a hole, 
is also visible in the saturability of the exciton-photon coupUng, resulting in an 
anharmonic interaction term which adds to the usual harmonic one: 
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dr[wi{r)Wc{r) + W*c{r)Wx{r)] 



1 



Wx[r) 

nsat 



(5) 



where Wsat = 7 /{16na^) is the exciton saturation density ||35]| . In GaAs, Wsat — 
2842 {jj.m)^^ and for a Rabi splitting of Qr = 4.4 meV, the ratio between satu- 
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ration and exciton-exciton interaction strength. 



^gxns, 



:0.1 



allow us to neglect the anharmonic term in (|5]l for the kind of physics we want to 
describe in this review. 

Therefore, the mean-field evolution of the coupled cavity photon-exciton dy- 
namics is described by the following non-linear Schrodinger equation or Gross- 
Pitaevskii equation (GPE): 



id, { "^^ 
\Wc 



Ho 



'gxlWxl^ 
Vcir] 



"^^ I . (6) 

Wc' 



Here, we have also added an external potential Vc{r) acting on the photon compo- 
nent, which later on we will use to describe the effect of photonic disorder present 
in the cavity mirrors. Note that Eq. (|6]) is a classical field description, which as- 
sumes the macroscopic occupation of a finite number of states, each described by a 
complex classical function y/. 



3 Optical parametric oscillator regime 

An accurate control of the polariton dynamics can be achieved by directly inject- 
ing polaritons at a given wave-vector and frequency with a properly tuned external 
laser — the resonant excitation scheme. In within this scheme two regimes can be 
singled out: (i) the regime where only the polariton state generated by the pump is 
a stable configuration of the system (we refer to this as the pump-only state); (ii) 
the regime where the polaritons continuously injected into the pump state undergo 
coherent stimulated scattering into the signal state (close to the normal direction) 
and the idler state (on the other side of the pump). Parametric scattering from pump 
to signal and idler can be self-induced by the continuous-wave (cw) laser above a 
pump strength threshold, in which case one refers to the optical parametric oscillator 
(OPO) regime. However, below the threshold for OPO, a second weak probe beam 
shined close to either the (expected) signal or idler states, can be used to 'seed' the 
parametric scattering processes and amplify the probe; in this case, one refers to the 
optical parametric amplification (OPA) regime. 

We introduce the concept of polariton parametric scattering and review the main 
experimental results on optical parametric amplification in the next section. In 



Sec. 3.2 we use a simplified theoretical model in terms of plane waves, for both 
the pump-only resonant state and the OPO state, summarising the main properties 
of both regimes and drawing an analogy with equilibrium weakly interacting Bose- 



Einstein condensates (BECs). Finally, in Sec. 3.3 we explain the necessity for carry- 
ing out a numerical analysis of the OPO. Much experimental work has been carried 
out on polaritons in the OPO regime lIl|36l[l9lll7l|38l|39l|40l|4T]|42l|43]|Bl|4a 
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(for a review on the experiments, see Ref. Il46l ). We will discuss the experimental 
achievements along with the theoretical description. 



3.1 Polariton parametric scattering and optical parametric 
amplification 

In the parametric scattering process, two polaritons from a pump mode, with wave- 
vector and frequency {kp, (Op}, scatter into a lower energy signal mode {kj, cOs} and 
a higher energy idler mode {k,, o,}. This scattering process has to conserve energy 
and momentum, therefore requiring that 

2kp = k, + k; 2cOp = (0, + (Oi . (7) 

This condition cannot be satisfied by any particle dispersion, for example paramet- 
ric scattering is forbidden for particles with a quadratic dispersion. In order to check 
whether parametric scattering is allowed for polaritons, one has to verify if the con- 
dition 

ICOLpikp) = 0)Lp{ks) + C0Lp{\2kp - k,|) (8) 

can be satisfied. If kj — 0, then the momenta of pump and idler are uniquely se- 
lected (see left panel of Fig. l2]l. In this case, the value of the pumping angle is also 
referred to as the "magic angle", and is located close to the inflection point of the 
LP dispersion. However, for a generic signal wave-vector, k; = k = {kx,ky), then for 
a fixed pump angle kp (assumed to be oriented along the x-direction, {kp,0), in the 
right panel of Fig. |2]i, the final states allowed in the parametric scattering process 
describe a figure-of-eight in momentum space B7ll28ll48l . 

In the case of optical parametric amplification experiments, parametric scat- 
tering is stimulated by a weak additional probe field. OPA was first observed 
in an InGaAs/GaAs/AlGaAs microcavity fWl, where a substantial signal gain 
of up to 70 was measured. Much experimental work has followed this first re- 
sult |l50l|5l]|52l|53l|5l|55l|56l|57l|58l|59l. Pump-probe parametric amplification 
of polaritons with an extraordinary gain up to 5000 and at temperatures up to 120K 
has been reached in GaAlAs-based microcavities and up to 220K in CdTe-based 
microcavities 155). In three-beam pulsed experiments lISTl . polaritons scatter from 
two equal and opposite angles, kp and — kp, into the LP and UP states at k = — 
note that at zero detuning, 5=0, 2co^ — COlp{0) + COupiO). Interestingly, parametric 
amplification has been also obtained for ultracold atom pairs confined in a moving 
one-dimensional optical lattice |60l. The role of the periodic optical lattice is to de- 
form the atom dispersion from the quadratic one, allowing parametric scattering to 
happen. 

The stimulated scattering regime can be reached also in the OPO configuration, 
i.e., without an additional probe beam. Now, stimulated scattering is self-initiated at 
pump powers above a threshold intensity, where the final state population is close 
to one. We will see that in this case, there is no special significance of the "magic 
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Fig. 2 Illustration of the basic idea of parametric scattering. Left: two LPs scatter from the pump 
state {kp, tt)p} towards the signal {kj , tt),} (here at zero momentum) and the idler state {k, = 2kp — 
kjjO); = 2 cQp — fflj} (at higher momentum), conserving momentum and energy. Right: Following 
Refs. (1211^, we plot |(»lf(^) + «L/'(|2k,; - k|) - 2fi)z.p(i(:,,)| as a function of k = {k^,ky). The 
white line is the zero value of the contour and the pump is oriented along the x-direction, k,, = 
(^p,0). The parameters used in both panels are the same as the ones of Fig. [T] 



angle", rather a broad range of pumping angles larger than a critical value, (0p > 10° 
for the parameters in Ref. ||6T| ) allow OPO with a signal emission close to normal 
incidence, 0j ~ 0°. In addition, for finite size pumping (see later Sec. 3.3 i, the pump. 



signal, and idler momenta are smeared in a broad interval, while their frequency still 
satisfy the matching conditions (|7| exactly. In the next three sections, we will focus 
mainly on the theoretical description of polariton resonant excitation with a cw laser 
field, describing the properties of first the pump-only state and then the OPO state. 



3.2 Bistability and OPO in the plane-wave approximation 



The theoretical description of polaritons in the resonant excitation regime can be 
formulated in terms of the same classical two-field non-linear Schrodinger equation 
previously introduced in Eq. (|6]l, where now an external driving field Fp{r,t) is 
added in order to describe the coherent injection of photons into the cavity: 



idt 










Ho + 



8xWx\ 







Vc{r) 



Wx 



A continuous-wave (cw) pumping laser can be written as 



(9) 



(10) 
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where -^fp.Op {r) can either describe a homogeneous pump with strength fp, -^f^.a,, {r) 



fp, or, as we will assume later in Sec. 3.3 a Gaussian or a top-hat spatial profile with 



strength fp and full width at half maximum (FWHM) Op. 

For a homogeneous pump, •^f„.a,,{r) = fp, and for a clean system, Vc{r) — 0, 
the conditions under which a stable OPO switches on can be found by making use 
of an analytical treatment ESl |6T1 l62l l63l . In fact, in this limit, each mode can be 
approximated as a plane wave. To simplify the analytical expressions, it is useful 
to rotate Eq. (|9| into the LP and UP basis, as described by Eq. Q, and to neglect 
the contribution from the UP states, assuming that the LP and UP branches are not 
mixed together by the non-linear terms. In this case, working in momentum space, 
VLp{r,t) ~ Lk^''"^V'z./',k(f)' Eq. (|9]) can be written as 

idt\\lLP^^ = {(JiLp{k) - iK{k)] \(/LP.k + 

L ^k,ki,k2 ¥lp.ki+k2-k¥LP.ki Wlpm + sin 0-t/pe"'™"'4,k^ , (1 1) 

ki.kj 

where K{k) = Kjf cos^ 0^ + Kc sin- 9^ is the effective LP decay rate and the interac- 
tion strength now reads as gk.ki.ki = gx cos Oi^ cos 9\]^^^]f^^_]^\ cosQ^^ cosQ^^^- 

If we consider solutions of the Eq. ( [TT| i where only the pump mode, k = kp, is 
populated, we can find an exact solution in the form 

xi/Lp{r,t) = y^e'C'pr-^.O y^^p^^^t) = pSkXpe-"""' , (12) 

where the complex amplitude p is given by 

[cOLp{kp) - (Op - iK{kp) +gx cos"* dkp IpP] P + sin Q^Jp = . (13) 

For practical purposes, one can substitute gxcos'^Oi^_^ i-^- 1 by redefining the pump 
strength /' = y^cos^ 9k sinOkfp, and rescaling the field strength p by \p'\ = 
y/gxcos^ 9k,\p\- Note that the ;if'^^' -non-linear interaction term, \p'\^p', renormalises 
the effective detuning of the pump mode from the LP dispersion, 

Ap = C0p-C0Lp{kp)-\p'\^ , (14) 

which now includes the blue-shift of the LP dispersion due to interactions. 

The general solution of the cubic equation ( [T3] l is well known, and exhibits a 
qualitatively different behaviour depending whether the pump frequency is blue- 
or red-detuned with respect to the LP dispersion. In particular, if cOp — C0Lp{kp) < 
V^K{kp), the system is in the optical limiter regime, where the population \p'\^ 
grows monotonically as a function of the pump intensity /'. If instead cOp — 
(OLp{kp) > \^K{kp), the system displays bistable behaviour, with a characteris- 
tic 5-shape of \p'\^ as a function of f' the second turning point coinciding with 
the point where the effective detuning Ap ( [T4] i changes sign (see Fig. l3]l. Because 
the branch with negative slope is unstable, the polariton density in the pump-only 
mode follows a hysteretic behaviour: Increasing the pump intensity, eventually the 
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Fig. 3 Dimensionless LP population |p'|^/K:(^p) = gxcos^'^i- |p|^/K:(^p) as a function of the 
dimensionless pump intensity /' = ^Jgxco'^ %_ sinO^ /,, for different values of the parameter 
5 = \(ap — aLp{kp)]/K{kp). When S < \/3 (left panel), the system is in the optical limiter regime, 
while for sufficiently blue-detuned pump frequencies, 5 > V^, a bistable behaviour is obtained 
(right). The sign of the interaction renormalised detuning, Ap l |14[ l, is also given. 



pump-only mode jumps abruptly into the upper branch, while if the intensity is then 
decreased, the polariton population decreases and jumps back down to the lower 
branch for smaller values of the pump intensity. Optical bistability in microcav- 
ity polaritons has been observed in Refs. fS?, "651, with evidence of a hysteresis 
cycle of the polariton emission as a function of the pump intensity. Multistability 
of two different polariton states, generated by either populating two different spin 
states Il56ll67l l68l or by injecting two states with two different pumps r69l has been 
also recently proposed and, in the spin case, observed. Part of the interest in this 
field is to realise all-optical switches llTOl and memories. 

The dynamical stability of the pump-only solution ( |T2] i can be established by 
allowing other states than the pump (i.e., the signal and idler states) to be perturba- 
tively populated via parametric scattering processes, 

VLPMit) = pSk^kpe-'""' + sSkMp-^e-''^""-''^' + /*4,k,+qe-'<™"+'*'*" , (15) 

where, {kv^,- = kp =f q, cOsj = Op =F Co}, and by assessing whether the time evolution 
of these states grows exponentially in time or not. Expanding to the first order in s 
and i, one obtains an eigenproblem for the amplitudes s and ; 



.2„2 



CO - A, - iK{ks) gxCsCiC^p 
gxcscicjp*^ -CO-Ai + iK{ki)J \i 



0, 



(16) 



where A,j = (O,, - (OLp{K^i) -2gxc]f'l\p\^ and Cp,.,,,- = cos0i,^„i,,i.. The complex 
eigenvalues O) can be obtained imposing that the determinant of the matrix in ( [T6] l 
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is zero. The dynamical stability is ensured if 3 (ft)) > 0. Therefore, the threshold 
for instability of the pump-only solution ( [T2] i can be found imposing the condi- 
tion 3 (ft)) = 0. By fixing the pump wave-vector and energy {kp,(Op) and the signal 
wave-vector k, (as well as the exciton and photon lifetimes, Kx.c), this provides a 
criterion for establishing the boundaries of the instability region, i.e., the lowest and 
highest values of the LP population \p\^ for which the pump-only solution is not 
stable. As shown in Refs. |61, 621^31, one can classify the instability as a single 
mode instability when q = and therefore kp = k, = k, — the Kerr instability. 
In particular, the branch with negative slope of the bistable curve (dashed line in 
Fig. [3| is always single mode unstable 116111631 . If instead q 7^ 0, the instability is 
parametric -like. Now, the total extent of the instability region corresponding to dif- 
ferent values of k^ is significantly larger than just the branch with negative slope. In 
addition, the OPO state does not require a bistable behaviour and can turn on also in 
the optical limiter case. In particular, it is possible to plot a "phase diagram" |[6TI of 
pump energy cOp as function of pump wave-vector kp, showing the regions where a 
pump-only solution is always stable, where the OPO switches on, and where instead 
a Kerr-type instability is only possible. In this way, in Ref. 1611 . it was shown that 
there is no particular significance to the "magic angle" for the pump. Rather, OPO 
conditions can be found for all angles larger than a critical value, 9p > Be (^ 10° for 
the parameters of Ref. lET]), as also confirmed experimentally Il39ll40l . In addition, 
the energy renormalisation of the polariton dispersion due to interactions moves the 
emission angles for the signal always close to 6^ ^ 16111401 . This is also confirmed 



by the numerical simulations we have carried out and illustrated later on in Sec. 3.3 
The method described above implies negligible populations of the signal s and 
the idler / and therefore it allows to find the conditions for the OPO threshold. In 
order to find the OPO states, one cannot linearise in s and /, but instead include 
the contributions of finite signal and idler populations to the dispersion renormal- 
ization (|61 1. In this way, in the region unstable for parametric scattering determined 
with the method described above, one can describe first the increase (switch-on) and 
later the decrease (switch-off) of the signal and idler populations as a function of 
the pump power. It is interesting to note that by doing that, i.e., by substituting ( fTS) 
into ( fTTj l, "satellite states" oscillating with energies ft)^^ — 2(0s ~ (Op — (Op — 2(0 and 
ftj/j = 2ft), — (Op = (Op + 2ft) automatically appear. In fact, above OPO threshold, 
when signal and idler populations are not negligible, parametric scattering from the 
signal (idler) state into the pump and second-signal (second-idler) satellite state take 
place, i.e., 2s i->-p + S2 i2i>-^ p + h) and therefore 2 ft)j == ftjp + ft)j, (2ftj, — (Op + at^)- 
This is clearly seen in the "exact" OPO solution obtained numerically (see, .e.g., 
Fig.|4]i, as well as it has been observed experimentally (see, e.g., Ref. Il37l ). One has 
to note however that the population of the "satellite states" by multiple scattering 
processes is always negligible w.r.t. that one of pump, signal, and idler (see right 
panel of Fig. HI. 

We will introduce the numerical modelling used to describe the problem for a 



finite-size pump later in Sec. 3.3 Before doing that, in the next section , we con- 
centrate on the analogies and differences between an OPO state and an equilibrium 
weakly interacting Bose-Einstein condensate (BEC). 
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3.2.1 Spontaneous t/(l) phase symmetry breaking and Goldstone mode 

The OPO state looks at first sight very different from an equilibrium weakly inter- 
acting BEC. In particular, the OPO is an intrinsically non-equilibrium state charac- 
terised by the (macroscopic) occupation of three polariton states only, one directly 
populated by the external pump and the signal and idler states populated by para- 
metric scattering. Contrast this with the thermodynamic phase transition in a BEC, 
where the macroscopic occupation of the ground state occurs when, for a thermal 
distribution of bosons, either the temperature is lowered below a critical value or the 
density is increased. The OPO state does however share with a BEC the fundamen- 
tal property of spontaneous symmetry breaking of the phase symmetry ISTl ITTI . In 
fact, the external laser fixes the phase of the pump state 0p and parametric scattering 
processes constraint the sum of the signal and the idler phase only, 2(j)p — (j)^ + (j)i, but 
leaves the system to arbitrarily choose the phase difference (j)^ — (j)i. In other words, 
one can easily show that the system of three equations one obtains by imposing the 
OPO solution ( [T5] l into the mean-field equation ( fTT) is invariant for a simultaneous 
phase rotation of both signal and idler states: 

s ^ se'^ i ^ ie-'^ . (17) 

This £7(1 ) phase rotation symmetry gets spontaneously broken in the OPO regime, 
where the signal and idler spontaneously select their phase, though not indepen- 
dently. Note that in this respect the OPO regime differs very much from the optical 
parametric amplification (OPA) regime, where both signal and idler phases are fixed 
by the external probe, and therefore the U{\) phase rotation symmetry is explicitly 
broken by the probe and no phase freedom is left in the system. 

Goldstone's theorem states that the spontaneous symmetry breaking of the L'^(l) 
phase symmetry in OPO is accompanied by the appearance of a gapless soft mode, 
i.e., a mode ft)(k) whose both frequency 9{[a)(k)] and decay rate 3[(a(k)] tend to 
zero in the long wave-length k — 7> limit. The dispersion for the Goldstone mode in 
OPO has been derived in Ref. ITTTI . where also an experimental set-up to probe it's 
dispersion has been proposed. In addition, the appearance of spontaneous coherence 
in OPO have been shown via quantum Monte Carlo simulations |72 1 through the di- 
vergence of the coherence length when the pump intensity approaches the threshold. 
In contrast, in the OPA regime, where the phase rotation symmetry is explicitly bro- 
ken by the probe, there is no Goldstone mode and a gap opens in the imaginary part 
of the elementary excitation dispersion, 3[G)(k)]. 

We would like to stress here that, even though an equilibrium weakly interacting 
BEC and an OPO state share the fundamental property of spontaneous symmetry 
breaking of the phase symmetry, some care needs to be applied in pushing this anal- 
ogy further. In particular, the existence of a free phase alone is not sufficient to 
ensure the paradigmatic properties of a superfluid, such as the Landau criterion, the 
stability of quantised vortices, and the persistency of metastable flow. For example, 
let us consider here the case of the Landau criterion: In an equilibrium weakly in- 
teracting BEC, the existence of the soft Goldstone mode (the Bogoliubov mode). 
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with it's characteristic linear dispersion for k — > 0, ft)(k) ~ c^k , implies the exis- 
tence of a critical velocity, Vc = mm]iCo{k) / k — Cs (the speed of sound), below which 
a perturbative defect dragged through the fluid cannot dissipate energy (superfluid 
regime). In the non-equilibrium OPO regime instead, similarly to what happens for 
incoherently pumped polaritons condensates Il73ll74ll2n . the unusual form of the 
excitation spectrum — diffusive at small momenta — poses fundamental questions 
on the fulfilling of the Landau criterion and the possibility of dissipationless super- 
flow. 

Similarly, properties such as the appearance and stability of quantised vortices 
and the persistency of metastable flow need to be independently assessed in po- 
lariton fluids in the three different pumping schemes available — (i) non-resonant 
pumping; (ii) parametric drive in the optical-parametric-oscillator regime; (iii) co- 
herent drive in the pump-only configuration. In fact, in the case of an equilibrium 
condensate, the ground state is flow-less, i.e. a vortex solution is unstable in non 
rotating condensates ^ . In contrast, in a polariton fluid, its intrinsic non-equilibrium 
nature implies the presence of a flow even when a steady state regime is reached. 
In this sense, not always the presence of vortices can be ascribed to the superfluid 
property of the system. We will discuss these aspects more in depth later in Sec. [5] 



3.3 Numerical modelling 



We have seen in Sec. 3.2 that, for homogeneous pumps, ^f.,.c„ = /p. the condi- 
tions under which a stable OPO switches on can be found analytically by assuming 
that pump, signal, and idler states can be described by plane wave fields ( [T5] l and 
therefore are characterised by single wave-vectors kp j,, and by uniform currents, the 
intensity and direction of which are given by k^.s , ^. However, for pumping lasers 
with a finite excitation spot, ^/„,(7„('")^ such as the ones employed in experiments, 
one can only resort to a numerical analysis f78 1 of the coupled equations (|9]). A fi- 
nite size pump implies that, in the OPO regime, pump, signal, and idler states are 
broaden in momentum; as a consequence, these states are going to be characterised 
by non-trivial configurations of the currents (see Fig. |5]l. We will see later on that 
these currents play a crucial role in the occurrence and dynamics of both metastable 
and spontaneous vortices in OPO. 



^ In rotating condensates, a vortex can be created if the angular velocity is higher than a critical 
value I75II76I . When rotation is halted, then the vortex will spiral out of the condensate 1771 . 
"^ Given a complex field or wave-function, \\l>{rj)\e'^^^''\ describing either a quantum particle of 
mass m or a macroscopic number of particles condensed in the same quantum state, the current is 
defined as 1761 : 

j(r,0 = -|v.(r,OpV<|.(r,0 = |v/(r,f)|'v.v(r,f) , (18) 

m 

following, with a slight abuse of notation, we will refer to the current as the gradient of the phase 

only, V0(r,f). 
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In particular, we numerically solve Eqs. (|9]l on a 2D grid of typically N x N = 
2^ X 2^ points and a separation of 0.47 jim (i.e., in a box LxL= 140 jUmx 140 ^m) 
by using a 5*-order adaptive-step Runge-Kutta algorithm. We have checked that our 
results are converged with respect to both the resolution in space L/N and the one 
in momentum n/L. Note also that of course the extension of the momentum box 
^max ~ 7tN/L has to be big enough to contain the idler state. In the specific case of 
Figs. [4] Is] andl6] we have chosen a smoothed top-hat profile •^fp.a,, ('") with FWHM 
Op = 7(J^m and (maximum) strength fp (later for Fig. [8] we have chosen instead a 
FWHM Op = 35 /im). Considering the case of zero detuning, 5 = 0, we pump at 
kp = 1.6 livcT^ in the x-direction, kp = [kp,0), and at cOp — 0)^ = —0.44 meV, i.e. 
roughly 0.5 meV above the bare LP dispersion, and gradually increase the pump 
strength until the OPO switches on. We find that broader LP linewidths imply a 
wider range in pump strength of stable OPO and for this reason we fix iCx = Kc = 
0.26 meV in these particular runs. We define /* as the pump strength threshold for 
OPO emission — here and in the following, we only select OPO solutions which 



reach a dynamically stable steady state (dynamical stability is studied in Sec. 3.5.1 1 
In the case of Figs. HI and [5] the pump strength is fixed just above threshold, /,, = 

1.25/J''- 
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Fig. 4 Left panel: OPO spectrum |v/c.x(k, (o)\^ for a top-hat pump of FWHM Op = 70 /im and 
intensity fp = 1.25/* above the threshold pump power for OPO, fp . For this particular nan we 
resonantly pump at kp = 1.6 /im^' in the x-direction, k,, = (kp,0), and at (Op — (O^ = —0.44 meV. 
Polaritons at the pump state undergo coherent stimulated scattering into the signal and idler states, 
which are blue-shifted with respect to the bare lower polariton (LP) dispersion (green dotted line) 
because of interactions. Cavity photon (C) and exciton (X) dispersions are plotted as gray dot- 
ted lines. Above threshold, as discussed in the text, we observe the population of the satellite 
states in addition to the one of signal and idler. Right panel: Momentum integrated spectrum, 
Ek IV'c,x(k, C0)|^, as a function of the rescaled energy ft) — ft)^. Pump, signal, idler and satellite 
states are all equally spaced in energy by roughly 1.19 meV. 
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The numerical analysis provides the time evolution of both photon and exciton 
fields either in space, y/c^xi^,t), or in momentum, v/cx(k,f )■ The OPO implies the 
simultaneous presence of (at least) three states emitting at different momenta, and 
therefore, at a fixed time t, the full emission v/c.x(r,0 is characterised by interfer- 
ence fringes. Because, as for the pump, the dominant wave-vectors for signal and 
idler are in the jc-direction, the fringes are vertical, i.e. predominantly oriented along 
the y-axis (see first panel of Fig. |5]l. We plot the photon component only, which 
is what can be measured experimentally. Note, however, that in cw experiments 
emission is always integrated in time, which clearly washes away the interference 
fringes. The OPO phase information can instead be recovered by obtaining inter- 
ference fringes with a reference beam in a Michaelson configuration. In addition 
to the full emission, either in space or momentum, one can also evaluate the spec- 
trum resolved in momentum \//cjf(k, Co) by taking the Fourier transform in time of 
Vc,x(k,f) (in FigWl a grid in time of 2^ points spaced by 0.3 ps has been used). As 
shown in Fig. H for the chosen parameters, a signal at 0)^ — O^ = — 1 .66 me V and an 
idler at a>s — (O^ = 0.75 meV appear with a sharp S-like emission in energy, which 
satisfies exactly the energy matching condition (|7]i, 2(0p = (Os + (Oi, as clearly shown 
by the momentum integrated spectrum on the right panel of Fig. [4] In contrast, the 
momentum distribution is broad (because of the pump being finite size) and peaked 
respectively at k^ ~ —0.2 /im^' and kj ~ 3.5 /im^', which only roughly satisfies 
the momentum matching condition, 2kp = ky +k,-. Note that the idler intensity is 
always weaker than the signal one because of the small photonic component at the 
idler. Further, note that, in addition to signal and idler, the spectrum also shows the 
appearance of satellite states {s2,sj,,. . . and (2, '3, ■ ■ ■ ) all equally spaced at around 



1.19 meV one from the other As discussed at the end of Sec. 3.2 their presence 
is a consequence of the secondary parametric scattering processes 2s i-^ S2 + p, 
2i H> /2 + p, 2*2 I— > ^2 + ss, and so on, which trigger on automatically as soon as 
signal and idler have finite populations. The occupation of the satellite states gets 
gradually suppressed the further we move higher in energy above the idler and lower 
in energy below the signal — which is why they are usually neglected in the plane 



wave approximation, as discussed in the end of Sec. 3.2 Note also that the satellite 



states just described do not imply the presence of phase symmetries additional to 



the t/(l) one described in Sec. 3.2.1 These satellite states therefore differ from the 
states which one could obtain as a result of secondary instabilities, e.g. 2s i-^ S2 + s'^ 
with S2 7^ p and 2/ h^ /j + i^ with i!^ ^ p, and successive spontaneous symmetry 
breaking mechanism |79|. 

In order to analyse the OPO properties, similarly to what is done in experiments, 
it is also useful to filter the full emission in order to select only the emission com- 
ing from the signal, pump or idler This can be equivalently done either filtering in 
momentum space in a cone around the momenta kp ^ ,■ or filtering in energy, bring- 
ing to the same results. We indicate the filtered spatial profiles of pump, signal, and 

idler by \\)f^'^^'{r,t)\e'^'^-'' ^'"' . The associated currents, V0^''^'', are a superposition 
of a dominant uniform flow kp j , (which is subtracted from the images of the sec- 
ond panel of Fig. Bll and more complex currents (caused by the system being finite 
size), which move particles from gain to loss dominated regions. Note that because 
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Fig. 5 Full emission (first panel) and filtered emission of signal (second), pump (third, with super- 
imposed cuiTents), and idler (fourth) states for the same parameters of Fig. [4] We plot the rescaled 
currents of the signal in the second panel by subtracting the dominant uniform flow, i.e., W(j)^ ^ — kj. 



we select only steady state OPO solutions, the profiles of pump, signal, and idler, 
|i/^'^'(r,f)|, are time independent. In addition, note that, the fact that the pump is 
shined on the microcavity with a finite angle respect to the normal incidence, im- 
plies that, for rotationally symmetric pump profiles, the symmetry inversion r i~-> — r 
is broken in the direction of the pump wave-vector k^. For example, if the pump is 
shined on the x-direction, kp = {kp,0), as in the case of Fig. E\ only the symme- 
try J '-^' —y is left intact. Clearly, this symmetry, while allowing vortex-antivortex 
pairs, does not in principle permit OPO solutions carrying single vortices, which 
can spontaneously appear in presence of a symmetry breaking perturbation, such as 



disorder (next paragraph) or a noise pulse (see Sec. 3.5. 1 1. 

The typical changes of the signal space profile as the pump power is increased 
above threshold, together with the pump, signal, and idler intensities, are shown in 
Fig. [6] For these runs we fix the parameters, such as 0)^, kp, and the pumping spot 
size dp, as in Fig. l4| but we also include a static photonic disorder potential — 
see Eq. (|6]l. In particular, here we consider a disorder potential with zero average, 
(Vc(r)) — and a spatial distribution. 



{Vc{r)Vc{r')) - a]e 



2 -|r-rf/2^2 



(19) 



with a correlation length 4/ — 20 fim and strength a^i ~ 0.1 meV. Below threshold 
the system is in a pump-only state. By increasing the pump power fp, above thresh- 
old, the OPO signal first switches on only in a small (compared with the pump 
spot FWHM C7p = 70 jUm) region (see inset 1). At fp = 1.2/* (inset 2) the signal 
becomes large and quite homogeneous, though, already at fp = 2.3/*, the OPO 
signal starts switching off in the middle (inset 4), and then it slowly switches off 
everywhere. A similar behaviour has been found in the numerical simulations of 
Ref. Ii78l (though there a small pump beam of FWHM Gp ^ 20 fim has been used), 
as well as observed experimentally in Ref. Il42l . 

The qualitative behaviour of the signal (as well as the idler) profiles, in particular 
their switching on and then off, as a function of the pump power that we have just 
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Fig. 6 Evolution of the signal, idler and pump state intensities as a function of the pump intensity 
fpl fp (Isft)- Space profiles of the filtered signal at different values of the pump intensity (right). 
The parameters are the same as in Fig. R] with the addition of a photonic disorder potential Vc (r) , 
correlation length lij ~ 20 /im and strength Oii ~ 0.1 meV (contour-level lines in the panels on the 
right). 



described for a disordered sample is very similar to the case of an OPO in a homo- 
geneous sample ^ (i.e., with no photonic disorder, Vc(r) = 0). One of the main dif- 
ferences is that for homogeneous samples the profiles are y H> — y symmetric, while 
this symmetry is explicitly broken by the photonic disorder In addition, the fun- 
damental difference between the homogeneous and the disordered case, is that the 
presence of photonic disorder promotes stable vortex solutions in large pump spot 
OPOs at intermediate pumping strengths, fp ~ 1 .4/* — such as the one shown in 
panel 3 of Fig.l6]which carries two vortices. Single or multiple vortex solutions are 
generally not allowed in the homogeneous case because of the y n- — y symmetry, 
which instead only allows pairs of vortex-antivortex solution y^^—y symmetric. In 
large pump spots, such as the one of Fig.|6j vortex-antivortex solutions in the clean 
case tend to be dynamically unstable, i.e. easily destabilised by a weak noise pulse. 



while, as analysed later in Sec. 3.5 spontaneous vortex solutions in homogeneous 



cavities can be stabilised by a small pump spot (see Fig. ^ confining the vortex in- 
side. Note finally that, spontaneous vortices in disordered cavities with a large pump 
spots are not pinned into minima of the disorder potential, rather, as analysed in the 



Sec. 5.2.2 are the OPO steady state currents in the signal to play an essential role in 
the stabilisation of vortices. 



3.4 Vortex phase and profile 



Before moving on to describe the occurrence of stable vortices in OPO, and, later, 
the onset and dynamics of metastable vortices, let us briefly remind the definition of 
a quantised vortex in an irrotational fluid. In general, a quantised vortex with charge 



^ Note also that we find that the value of the pump threshold for OPO is not altered by the presence 
of a weak photonic disorder. 
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m is described by a wave-function, 



V/(r) ^ \(fo{r)e 



im(p[r) 



(20) 



the phase of which, m(p{r), Hnearly winds around the vortex core from to 2nm 
(with m integer) — i.e., in cyhndrical coordinates centered at the vortex core, (p is 
the azimuthal angle. This impHes that the vortex carries a quantised angular mo- 
mentum, hm. In addition, the phase has a branch-cut and therefore is not defined at 
the vortex core, implying the vortex wave-function has to be zero at the vortex core. 
An example of an m = — 1 vortex, with i/o('") = re^'"' '^^i ', has been plotted on the 
left panel of Fig.lT] Here, the phase winds clock- wise around the core, from to 2n, 
and therefore the vortex current. 



V(p(r) = 



(P 



(21) 



is constant at fixed distances from the vortex core, r, while decreases inversely pro- 
portional to the distance (right panel of Fig. |7]i. Contrast this with the case of a 
rotational vortex in a classical fluid which rotates as a solid body with an angular 
velocity Q: Now, the fluid tangential velocity is zero at the vortex core and increases 
linearly with the distance, i.e., \,p — Qr(p. Quantised vortices can be detected in in- 
terference fringes (middle panel of Fig. U) as fork-like dislocations, the difference 
in arms giving the charge \m\ of the vortex. 
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Fig. 7 Typical profile (left), phase and currents J2TJ (right), and interference fringes (middle) of 
an m = — 1 vortex <20k. 



3.5 Stable vortices in a small sized OPO 



As explained later in SecB] spontaneous stable vortices differ from metastable vor- 
tices (described in Sec. |5Tl[): Metastable vortices can only be injected externally. 
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e.g. by an additional Laguerre-Gauss beam probe, into an otherwise stable symmet- 
ric state, and their persistence is due to the OPO superfluid properties ||8]|22l- The 
metastable vortex is a possible but not unique stable configuration of the system. In 
contrast, as for non-resonantly pumped polaritons Il80l l5l. the appearance of sponta- 
neous vortices is not a consequence of the polariton condensate being superfluid, but 
rather to the presence of currents related to the non-equilibrium nature of these con- 
densates. This strongly differs from the case of equilibrium superfluids, the ground 
state of which is flow-less. Later, in Sec.lTlwe will briefly discuss how, for polaritons 
non-resonantly injected into a microcavity, the presence of a confining potential can 
generate currents favourable to the spontaneous formation of vortices |5, 81 1 and 
vortex lattices ifSOJ . 

For resonant excitation, currents arise in the OPO regime due to the simulta- 
neous presence of pump, signal, and idler emitting at different momenta, as well 
as by the fact the system is finite size (see Fig. Hli. We have seen in Fig. l6]that, 
similarly to non-resonantly pumped polaritons, the presence of a disorder potential 
can lead to the spontaneous appearance of vortices. However, it is remarkable that, 
even in the absence of disorder or trapping potentials, the OPO system can undergo 
spontaneous breaking of the y i— > — y symmetry and become unstable towards the 
formation of a quantised vortex state with charge m = ±1 if the size of the OPO is 
small enough |22|. This is the subject of this section. Further, as discussed in some 



detail later in Sec. 3.5.2 like for equilibrium superfluids, both stable and metastable 
vortices are characterised by a healing length which is determined by the parameters 
of the OPO system alone. Spontaneous stable vortex solutions are robust to noise 



(Sec. 3.5.1 1 and to any other external perturbation, and thus should be experimen- 
tally observable. However, while spontaneous vortex solutions in OPO have been 
observed for a toroidal pump spot ^, so far they have not been observed in OPO 
with a 'simply connected' pump profile, e.g., either a Gaussian or a top-hat. 

3.5.1 Dynamical stability 



As mentioned in Sec. 



3.3 



if the pump is shined on the x-direction, kp = {kp,0), only 
the symmetry y i— ?> —y is left intact in the system. Clearly, this symmetry, allows for 
OPO solutions where the signal (and therefore also the idler) have vortex- antivortex 
pairs, with the vortex core position at {xc,yc) and the antivortex core position at 
{xc,—yc)- However, both single and multiple vortex solutions explicitly break the 
y I—)- —y symmetry and cannot be accessed by the dynamics — Note that two vortices 
located at opposite sides with respect to the x-axis break the y i— ?► —y symmetry 
because of the currents. 

In order to check the dynamical stability of OPO states, one has to add small 
fluctuations to the steady state mean-field solution: The existence of modes with 
positive imaginary part in the excitation spectrum indicate dynamical instability 
towards the growth of different modes. The dynamical stability analysis for OPO 



D. Sarkar (University of Slieffield), private communication. 
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t=432ps t=1.2ns 




Fig. 8 Appearance of a spontaneous stable vortex solution in a homogeneous small sized OPO. 
Filtered signal profile |i|/^(r,?)| with superimposed currents V<p^{rj) (upper panels (a-c)) and full 
momentum emission \ij/c{kx,0,t)\ (lower panels (d-f), in arb. units) at three different times: f = 
(a,d), t = 432ps (b,e), and 1.2ns (c,f). At f = a pulsed weak random noise of strength 0.01 (see 
text) is added to the OPO steady state (in (d) both OPO momentum profiles without and with the 
added noise are shown for comparison) and at f = 432ps a vortex, with m = — 1 , enters the signal 
and settles into a steady state. Note that, because of phase matching conditions, the presence of 
an m = — 1 vortex in a signal implies the presence of an m = 1 antivortex in the idler. A vortex 
(antivortex) in the signal (idler) space emission appears also as a dip in momentum space at the 
signal (idler) momentum (e,f). Parameters used: smoothed top-hat pump with FWHM CTp = 35/im, 



pump strength fp = 1.12/, 



.(th) 



1.6 /im in the ;i:-direction, fi)„ 



-0.44 meV, zero 



detuning 5 = 0, and Kx = Kc = 0.22 meV. Adapted from L22i . 



described in within the plane-wave approximation of Sec. 3.2 has been discussed in 



Refs. IISTIItTI . Equivalently, stability can be numerically checked by introducing a 
weak noise. In particular, we add white noise as a quick (5-like in time) pulse at a 
certain time to to both modulus and phase of excitonic and photonic wavefunctions 
in momentum space, \\ifx.ciKt)\e"^'^''^^^''^'- 

lrx,c(k,fo)| 1-^ |v/x,c(k,fo)|+5|v/x,c(k)| 
fe,c(k,fo) i-^ |0x,c(k,fo) + 5(/)x.c(k) . 

Both 5|v/x,c(k)| and 50x,c(k) are white noise functions, with an amplitude In for 
the the phase 5(^x,c(k), while the amplitude of the noise in the modulus 5| V/x.c(k)| 
is specified in units of the maximum of the pump intensity in momentum space. 
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Following this procedure, we have been able to single out symmetric OPO states, 
as shown in Fig. ISla), which are unstable towards the spontaneous formation of 
stable vortex solutions. After the y >-^ —y symmetry is broken by the noise pulse, 
we have observed a vortex with quantised charge m — ±1 (m = =Fl) entering and 
stabilising into the OPO signal (idler) — Note that parametric scattering constrains 
the phases of pump, signal, and idler by 2(j)p = 0.s + 0,- (see Sec, |3.2.l"] i, therefore an 
m = —\ vortex in the signal at a given position implies an m — 1 antivortex in the 
idler at the same position and vice versa. In the case of Fig. [8] and the right panel 
of Fig. l9] the noise strength is 0.01 and 432ps after the noise pulse, a vortex with 
m = —\ (m = +l) enters the signal (idler) and stabilises. The strength of the noise 
has no relevance on the final steady state, and in particular it can be infinitesimally 
weak. Different noise strengths do only affect the transient time the system needs 
to accommodate the vortex and reach the new steady configuration. We have in 
addition examined whether the vortex steady state is dynamically stable by applying 
an additional noise pulse. For weak noise, with a strength up to 0.1, the vortex is 
stable and can only drift around a little before settling again into the same state. For 
strong noise, with strength 1 and above, the vortex gets washed away, but after a 
transient period, the very same state enters and stabilises again into the signal, with 
the possibility of flipping vorticity ^. Different noise strengths do not affect the final 
steady state, but only the transient time. 
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Fig. 9 Steady state filtered signal profile (dotted line) y/^(x,0,/) for y ~ before the airival of 
either a Laguerre-Gauss vortex probe \2A\ with a pi, ~ 1/im (left panel, red dashed line) or a noise 
pulse of strength 0.01 (right panel, red dashed line) — same OPO conditions as Fig. [8] After the ar- 
rival of any perturbation breaking the y^-^ —y symmetry, the same vortex with charge m = ± 1 (solid 
shaded curve) stabilises into the signal. [From [221 ask for copyright permission ! ] 



' When generated by a noise pulse, both stable and metastable vortices have equal probability to 
have either charge ±1. Similarly, when vortices are triggered via a Laguerre-Gauss probe, their 
vorticity can flip during the transient period. In particular, flipping can follow the appearance of 
two antivortices at the edge of the signal, one recombining with the triggered vortex. Note that 
the vorticity flipping conserves the total orbital angular momentum, in the sense that when for the 
signal m flips, say, froin +1 to —1, for the idler the opposite happens, i.e. m flips from —1 to -|-1. 
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As discussed later in Sec. 5] one can alternatively break the y i-> —y symmetry 
by a pulsed vortex probe (|24 1, and assess whether the stable steady state is in any 
way dependent on the external perturbation. The homogeneous OPO states which 
are unstable towards the spontaneous formation of stable vortices following a white 
noise pulse, exhibit the same instability following a vortex Laguerre-Gauss (LG) 
probe pulse (see the left panel of Fig. |9]l. The steady state vortex is independent on 
both the probe intensity fpt, and size Gpi,, however the weaker the probe the longer 
the vortex takes to stabilise, between 30 and 400ps for our system parameters. As 
shown in Fig. l9] the stable vortex following the LG probe is exactly the same as 
the one triggered by a weak white noise, indicating that the probe acts only as a 
symmetry breaking perturbation. 

Summarising, one can find OPO conditions where the y i~> — y symmetric so- 
lution is dynamically unstable and any symmetry breaking perturbation allows the 
signal and idler to relax into a stable steady state carrying a vortex with charge ± 1 . 
For homogeneous cavities, i.e., in absence of any disorder or confining potential, 
we found that this requires either a small Gaussian or small top-hat like pump spot 
which can confine the vortex inside or a doughnut-shape pump spot. Instability of 
the uniform state to spontaneous pattern (e.g., vortex) formation is a typical fea- 
ture of systems driven away from equilibrium |79 1. Similarly we find conditions for 
which the uniform OPO solution is unstable to spontaneous formation of a quan- 
tised vortex. In alternative, a disorder potential breaks the symmetry explicitly and 
allows the pinning of stable vortex solutions in OPO, which is less surprising. 



3.5.2 Healing length 

In contrast to their classical counterpart, quantised vortices with the same angular 
momentum \m\ are all identical, with a size (or healing length) determined by the 
system non-linear properties [76 1. In the case of a superfluid in equilibrium with 
a typical interaction energy gn (n is the average density) and mass m, the heal- 
ing length, ^ = l/y/2mgn, is the typical distance over which the condensate wave- 
function recovers its 'bulk' value around a perturbation. In particular, for an \m\ = 1 
vortex ( |20] l, ^ is the typical size of the vortex. 

Similarly, in OPO, one case show that, like in equilibrium superfluids, both sta- 
ble (see Sec. 3.5 i and metastable (see SecB} vortices are characterised by a healing 



length which is determined by the parameters of the OPO system alone. In partic- 
ular, shape and size of the metastable vortices described in Sec. [5] are independent 
on the external probe. In the case of vortices in OPO, an approximate analytical 
expression for the vortex healing length can be derived for homogeneous pump- 
ing II22I I9I. assuming that only signal and idler can carry angular momentum with 
opposite sign, ±m, i/A*''(r) = .Vn^Je'*'" '"e*""'''f *''(r), while the pump remains in 
a plane-wave state, il/^ir) — ^/n^e'^i''^, as also supported by our numerical analy- 
sis. For pump powers close to OPO threshold, it can be shown |22 9| that signal 
and idler steady state spatial profiles are locked together and satisfy the following 
complex GP equation 
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where |a| ~ gx\/njH- This equation describes a vortex profile 1761 with a healing 
length given by: 

t, = {ImcgxVrhmy^'^ ■ (22) 

This expression is similar to the one of an equilibrium superfluid, with the conden- 
sate density replaced by the geometric average of signal and idler densities. Further 
above threshold, one can show that signal and idler profiles are no longer locked 
together, and that they start to develop different radii. In both the simulations of 
Figs. M and [m we find ^ ~ 4/im, compatible with the estimate ( p2j l. 

In Ref. 19|, vortices in OPO have been created in a controlled manner by adding 
a weak continuous probe in resonance to the signal. Even if the phase freedom of 
the OPO system is explicitly broken in this configuration by the vortex cw probe, 
because the ratio of the probe to signal power density is low, the size of the vortex 
has been demonstrated to be determined by the OPO non-linear properties only 
rather than by the imprinting probe. In particular, a systematic study of the decrease 
of the vortex core radius with increasing pump power above threshold has allowed 
to confirm the behaviour described by the Eq. (|22]). 



4 Triggered optical parametric oscillator regime 

Before moving on into the description of metastable vortices in OPO, i.e. vortices 
which are transferred by a pulsed vortex probe into the OPO signal and idler, and 
their relation to superfluidity (Sec. |5]l, we describe here first the effect of an addi- 
tional pulsed probe on OPO in general terms. As described previously, in the OPO 
regime, polaritons are continuously injected into the pump state, and undergo coher- 
ent stimulated scattering into the signal and idler states. The OPO is a steady state 
regime, where the filtered profiles of signal, idler, and pump, \y/^"l^'{r,t)\, are time 
independent. This also is reflected in the typical flat dispersion around pump, signal, 
and idler which can be observed in the OPO spectra (see Fig. [4]). The group velocity 
of pump, signal, and idler, defined as the derivative of the energy dispersion at kp ^ ,, 
is therefore zero. This however does not mean that there is no flow of polaritons, 
which instead is described by the phase velocity or current, V0^'^'' (see footnote 2 



on page 13 i, with a dominant uniform flow given approximatively by kp 

In resonantly pumped polaritons, in order to initiate a travelling wave-packet 
characterised by a finite group velocity, one needs to use an additional pulsed laser 
beam on top of the cw one. The description of the system is therefore still in terms 
of the Eqs. (|9]l, with a total pump term given by the sum of the cw laser ( [TOJ l and a 
probe beam Fpi, (r, f ) : 

F{r,t)^Fi,{r,t)+Fphir,t). (23) 
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For the moment being we will consider the generic case of a pulsed probe with a 
Gaussian space profile, shined at a momentum and energy {kph, (Opt} ^'■ 

Fpbir, t) = /p,,e-|'-VlV(2<^,?„)g'(k„,-r-ffl,,r)^-(r-,,4)V(2<T2) (24) 

A pulse duration of 3 ps (defined as the FWHM in time of Fpi,{r,t)) corresponds to 
oy = 1.3 ps. The idea, first introduced by Ref. [201, is that the pulsed probe triggers 
parametric scattering ^ between the probe state at momentum and energy {kpi,, COpi,} 
and a conjugate state at {k^. = 2kp — k^/,, cOc = 2(0p — (Opb} — because one can 
either have k^^ > kp or kph < kp, we use the state labels 'probe' and 'conjugate', 
rather than 'signal' and 'idler'; by doing so, one also doesn't confuse the states 
generated by the OPO with the additional ones generated by the probe. Both probe 
and conjugate states are travelling decaying states which can evolve freely from the 
laser probe constraints once the pulse switches off. Such states are referred to as 
triggered-OPO (TOPO) states. Note that a TOPO can be triggered in two regimes: 
either (i) in a regime where the cw laser drives the system above threshold for OPO, 
in which case the probe and conjugate states are the extra population states on top 
of the steady-state OPO signal and idler states, or (ii) when no OPO is present, i.e. 
for the cw pump strength below threshold. For simplicity, the numerical analysis 



discussed below in Sec. 4. 1 is conducted in the regime (ii), but we have checked 
that the qualitative results also hold in the regime (i) — where, now, the steady 
state OPO population needs to be subtracted so that one studies the properties of the 
population triggered by the probe only. 



4.1 Theoretical description of the TOPO 

cross-refer to Fabrice Laussy' chapter In order to analyse the dy- 
namical evolution of a TOPO wave-packet, we study numerically the time-dependent 
solutions of the equations (|9]l, with a total pump given by ( |23| ) and ( |24] i. The probe 
triggers parametric scattering between a probe state and a conjugate state. In the 
majority of cases, as discussed in Ref. [23 1, the parametric scattering is too weak 
to induce any significant amplification, and an exponential decay of both probe and 
conjugate populations is observed immediately after the probe Fpi,{r,t) switches 
off Here, the spectrum shows a strong emission from the pump state and a weak 
emission from the LP states mainly at momenta kp^ and k^. 

However, we have found conditions under which both signal and conjugate states 
get initially strongly amplified by the parametric scattering from the pump, then de- 



^ Note that, differently from the cw laser beam, the energy distribution spectram of which is 
essentially 5-like, a pulsed beam has an intrinsic width in energy, proportional to the inverse pulse 
duration, c, '. 

' If the cw pump drives the system into the OPO regime, then the parametric scattering triggered 
by the pulsed probe will be in addition to the one related to OPO. However, as discussed later, the 
TOPO regime can be reached also in absence of the OPO. 
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Fig. 10 Spectrum (left) and spatial profiles of pump and filtered signal |i//^(x, 0, f)| (right) for the 
TOPO regime. A short, oy = Ips, m = 2 Laguerre-Gauss ([25j (left) or Gaussian m = J24j (right) 
probe shined at kp^ = (1.4,0) ^um^' triggers the propagating probe and conjugate states, which 
lock to the same group velocity (for these simulations we fix Kx = and Kc = 0.02meV). Adapted 
froml23l. 
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Fig. 11 TOPO signal intensity (black line) and the intensity of an external Gaussian probe (red 
line) as a function of time. Parameters are the same ones of the right panel in Fig. [10] The TOPO 
signal gets initially strongly amplified, then decays slowly, and finally exponentially. 



cay slowly and, only at later times, decay exponentially (see Figs. [TO]and[TT]l — 
we refer to this as the 'proper' TOPO regime. A similar behaviour has been also 
observed in experiments ^. Now, the spectrum is observed to be linear, (O = \g k 



(see Fig. 10 1. A linear spectrum can be explained by the fact that, in order to have 
efficient parametric scattering, probe and conjugate state must have a large spatial 
overlap and therefore similar group velocities. Thus signal and conjugate group ve- 
locities need to lock, which is only possible if the dispersion becomes linear — a 
similar result has been also found in ID simulations (see Fabrice Laussy' s 
Chapter), as well as in experiments ll20l . The group velocity is defined as the 



^ See, for example. Fig. 3 of Ref. |82|, where the intensity maximum of the extra population is 
reached within 4 ps after the maximuin of the pulsed probe, is followed by a slow decay. 
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derivative of the energy dispersion with respect to the momentum. However, we can 
also measure it from the probe and the conjugate density variations in time, i.e.. 



>,/^^.'' - ^^P''-'^ . 



~ph,c 



as Vg ' = dr,„ ' /dt, where r,„ ' is the maximum of either the probe or conjugate 
spatial profile, which we use as a reference. By analysing the change in time of 
the spatial profile of the TOPO probe state \\j/^'^{r,t)\, it is possible to show 

that its group velocity Vg is given exactly by the derivative of the lower polari- 
ton (LP) dispersion evaluated at kp/,, i.e., for zero detuning and low densities by 



v[^ = ki,b/i2mc)-ky{2mcJk^+4,nlQl) (see Fig 



Kpb 



sistent with the form of the spectrum shown in Fig. 10 
to determine 



12 1. This behaviour is con- 



Further, we have been able 



that the TOPO linear dispersion is tangential to the LP branch at 



kph, thus its slope is given in this case also by 
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Fig. 12 Left panel: Probe (solid lines) and conjugate (dashed lines) density profiles at different 
times after the arrival of the probe (same parameters as in Fig.|10[l. Right panel: Group velocity, 
Vg, of the propagating probe state as a function of the probe momentum kph. The black dots are 
determined from simulations, whereas the solid (red) line is the derivative of the LP dispersion 
evaluated at kpi,, I'f . The blue dashed line is a guide for eye to indicate where the LP dispersion 
deviates from the quadratic. Adapted from 1231 . 



From the PL spectrum we can also deduce the nature of the wave-packet propa- 
gation. For systems characterised by a linear dispersion, like in the TOPO regime, 
one expects a soliton-like behaviour, where probe and conjugate states propagate 
without changing neither their shape nor intensity. For quadratic dispersion, a Gaus- 



pb 



sian wave-packet moves at a constant velocity vj 

shape in time but its width grows (FWHM=((Jp^ + ( 2,„ ^ ^ 

ever, that, due to the finite polariton Ufetime, the total density decays exponentially. 



2mc 
t 



and it preserves its overall 
)2)i/2) 1761. Note, how- 



' In the regime where the probe generates only a weak parametric scattering, aside the strong 
emission from the pump state, the dispersion is simply that of the LP, and thus is not surprising 
that the signal propagates with a group velocity given by v'jf ■ Remember that here the cw pump is 
below threshold for OPO. 
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with a rate given by '^2'^"' ^^ ^^^'^ detuning. Finally, for non-quadratic dispersion, 
propagation becomes complex: The wave-packet gets distorted and there are beat- 
ings in the spatial profiles. In general, due to the dynamical nature of the TOPO 
state, the system evolves between these different scenarios. In particular, only in the 
strong amplification regime the spectrum is linear, while it evolves back to the LP 
one at longer times. However, for a one-dimensional version of the equations (|9]), 
and for uniform, infinitely extended in space, pumping spots, non-decaying, soli- 
ton solutions have been recently found ll83l . This has been also generalised to a 
two-dimensional infinite systems fS?!, where, for some narrow range of pumping 
strengths, a soliton-like behaviour has been predicted for kpt, = and kpi, = kp. 
However, to date, a non-decaying wave-packet propagation has not been found in 
experiments. 

Finally we would Uke to note that, the typical behaviour of probe and conjugate 



states (left panel of Fig. 12 1 is analogous to the one discussed in four-wave-mixing 
experiments |85, 86 1: when the probe arrives, and shortly after that, the conjugate 
propagates faster then the probe, before getting locked to it with a small spatial 
shift of their maximum intensities. At later times, when the density drops and the 
parametric process becomes inefficient, the two wave-packets start unlocking — the 
conjugate slows down with respect to the signal if kc < kpt, as in Fig. 12 or it moves 
faster when kc > kpt,. 



4.2 Experiments 

The TOPO regime has been recently studied in experiments in Refs. Il20ll45]| (for a 
review see Ref |87 1). As previously described, the additional pulsed probe has been 
used to create a travelling, long-living, coherent polaritons signal, continuously fed 
by the OPO. A large increase of the signal lifetime has been observed for a pump in- 
tensity approaching and exceeding the OPO threshold |45|. This observation can be 
explained in terms of a critical slowing down of the dynamics following appearance 
of a soft Goldstone mode in the spectrum close to threshold. It is also consistent 
with the nature of wave-packet propagation in systems with linear dispersion. This 
has been used to interpret subsequent experiments, where the linearisation of dis- 
persion leads to the suppression of weak scattering and therefore to a polariton mo- 
tion without any dissipation ||20l[82l cross refer to Fabrice Laussy' 
chapter. Due to the finite size of the excitation spot, the travelling TOPO signal 
lives only as long as it reaches the edge of the excitation spot. However, as discussed 
in detail in Ref |87|, in order to asses the sustainability in time of the TOPO process, 
both pump and probe beams can be chosen so that the probe state forms at kpi, 0. In 
such a case, the polaritons in the probe state are not travelling and therefore and it is 
thus possible to measure the lifetime of the TOPO wave-packet, which is of the or- 
der of a nanosecond. The decay of the TOPO population in time, as well as the finite 
lifetime of TOPO pulses, indicate that the soliton behaviour predicted in Ref. Il84l 
is not the explanation of the current experiments. However, the linearisation of the 
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system's dispersion due to the parametric process, as well as the appearance of the 
Goldstone mode, provide a sufficient explanation of dissipationless propagation in 
free space, as well as frictionless flow against an obstacle, during the part of the 
dynamics when parametric processes are strong and the spectrum linear. 



5 Triggered metastable vortices as a diagnostic of the OPO 
superfluid properties 

OPO condensates, as well as polariton condensates pumped incoherently, share with 
weakly interacting Bose-Einstein condensates at equilibrium phenomena like the 
spontaneous breaking of the phase symmetry and the appearance of a Goldstone 



mode (see Sec. 3.2.1 1. However, being intrinsically non-equilibrium, all polaritonic 
systems need continuous pumping to balance the fast decay and maintain a steady 
state regime. In strong contrast with equilibrium superfluids, the ground state of 
which is flow-less, pump and decay lead to currents that carry polaritons from gain 
to loss dominated regions. This can lead to the spontaneous formation of vortices: 
The presence of currents in polariton condensates can lead to the spontaneous ap- 
pearance of vortices without invoking any superfluid properties. This is true for 
incoherently pumped polaritons in presence of a confining potential ||5][80][8T1, as 
well as for polaritons in the OPO regime, with the difference that here, even in the 
absence of disorder or a trapping potential, the system becomes unstable towards the 
formation of a quantised vortex state with charge m — zLl (see Sec. |3.5| l. In addition, 
the hydrodynamic nucleation of quantised vortices can appear as a consequence of 
the collisions of a moving polariton fluid with an obstacle, as will be briefly dis- 
cussed in Sec. IT] Therefore, in general, for polaritonic systems, one has to apply 
some care when using the appearance of vortices as a diagnostic for the superfluid 
properties of such a non-equilibrium system. 

In the case of equilibrium superfluids, the rotation of a condensate is accompa- 
nied, above a critical angular velocity |76 , 75 1, by the creation of quantised vortices. 
Here, vortices are stable as far as the system is kept rotating and become unstable 
when the imposed rotation is halted f77l. However, persistent flow can be observed 
when a BEC is confined into a toroidal trap and the quantised rotation is initiated 
by a pulsed Laguerre-Gauss beam |88] [89] l90l . The toroidal trap is essential to 
allow the vortex stability, because of the energy cost of the vortex core to move 
through the high density region from the center of the torus where the density is 
zero. The very same idea of questioning the persistency of flow in a BEC via a 
pulsed Laguerre-Gauss beam as a diagnostic for superfluidity, can be applied to po- 
laritons '". As recently proposed for non-resonantly pumped polariton condensates 
in Ref. ll2Tl . this definition of superfluidity as metastable flow is equally meaning- 
ful in non-equilibrium systems as in equilibrium ones. However, as we will see, the 

'" Note, however, that even if in the atomic and polaritonic case the same Laguerre-Gauss laser 
field is used, the mechanism of spinning the BEC atoms is different from the one which rotates 
polaritons. 
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important difference is that, in the OPO regime, flow persistency is possible even in 
a simply connected geometry, i.e., without the need of a toroidal trap which pins the 
vortex. 

A pulsed Laguerre-Gauss (LG) probe beam carrying a vortex of charge m can be 
described by: 



^pi,(r,f)-/pfc|r-r,fc|l"V™'Pe- 



i'n<P„-\'^-^pb\^H^(ylb)j(\ip,,-r-(i>p,,t)„-{t-tpt?l(2a}) 



(25) 



with {kp/j, (Opt] can be tuned resonantly to either the OPO signal or idler As dis- 
cussed in the next section, by using a pulsed LG beam ( p5] l, vorticity has been shown 
to persist not only in absence of the rotating drive, but also longer than the gain in- 
duced by the probe, and therefore to be transferred to the OPO signal, demonstrating 
metastability of quantised vortices and persistence of currents in OPO. Experiments 
and theory will be discussed in the next Sec. 5. 1 




Fig. 13 Time evolution of the polariton signal following the arrival of a LG pulsed beam carrying 
a vortex of m = 1 . The first raw are the interference images obtained by overlapping the vortex 
with a small expanded region of the same image far from the vortex core, where the phase is 
constant, while the second raw are the space profiles of the signal. The sequence demonstrates 
that the vortex remains steady as a persisting metastable state for times much longer than the extra 
population created by the probe pulse and eventually gets imprinted in the steady state of the OPO 
signal. This is revealed by the strong contrast of the fork in the interference images for as long as 
the core remains within the condensate area. From Ref . L8J . 



5.1 Theory and experiments 



In the case of metastable vortex solutions, the symmetric vortex-less OPO steady 
state is dynamically stable, but, because of its superfluid properties, can support 
persistent metastable currents injected externally. From a theoretical point of view. 
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metastable solutions can be equally induced by either a vortex probe pulse ( |25| l or a 
noise pulse. However, differently from the case of stable vortices, metastable solu- 
tions require a threshold in the perturbation breaking the system y h-s- —y symmetry. 
For the simulations of Ref |8|, as we were interested mainly into the transfer of 
angular momentum from the probe into the OPO signal and idler, we have been 
considering conditions where the parametric scattering induced by the probe is too 
weak to induce any significant long-lasting amplification, and the gain introduced 
by the probe on top of the OPO disappears quite quickly. We have found conditions 
where the vortex is transferred from the probe into the signal ' ' (and antivortex in 
the idler) immediately '^ when the probe is shined. The transfer is followed by a 
transient time during which the imprinted vortex drifts around inside the signal and 
in certain cases settles into a metastable solution. Similarly to what happens to sta- 
ble vortex solutions, we have found that the spatial position of the metastable steady 
state vortices is close to the position where the OPO signal has the currents pointing 
inwards (see second panel of Fig. El. The influence of currents on the formation of 
vortices is discussed further in Sec. |5.2| Such metastable solutions do not always 
exist: if the probe is positioned well inside a wide OPO signal, as the creation of 



a vortex is accompanied by the creation of an antivortex (see Sec. 5.2 1, often, the 
vortex-antivortex pair quickly recombines; in other cases, during the transient pe- 
riod, the excited vortex can spiral out of the signal. Finally note that, as discussed 



in Sec. 3.5.2 the shape and the size of metastable vortices are independent on the 



external probe but are only determined by the parameters of the OPO. 



In the experiment of Ref. 1 8 1 also shown in Fig. 13a vortex is excited by a probe 
smaller than that of the signal to allow free motion of the vortex within the conden- 
sate. Vortices are detected, and their evolution in time followed by a streak camera, 
in interference images, generated by making interfere the OPO signal with a con- 



stant phase reference beam in a Michelson interferometer (second row of Fig. 13 1. 
As single shot measurements would give a too low signal to noise ratio, every picture 
is the result of an average over many pulsed experiments taken always for the same 
OPO conditions. The probe triggers a TOPO response, creating a strong gain and an 
extra decaying population on top of the OPO signal (TOPO). In experiments, dif- 
ferent regimes have been investigated. In particular, it has been possible to establish 
that, only under very high pump power and at specific points in the sample, the vor- 
ticity was transferred from the TOPO into the OPO signal, generating a metastable 
vortex solution. This not only demonstrates that the OPO polariton condensate can 
show unperturbed rotation, but also that a vortex can be another metastable solu- 
tion of the final steady state, demonstrating therefore the superfluid behaviour in the 
non-equilibrium polariton OPO system. After the vortex is imprinted into the OPO 



^' We checked that m = ± 1 (m = ^ 1) vortex solutions can appear only into the OPO signal (idler). 

A vortex probe pulse of any charge m injected resonantly to the pump momentum and energy gets 

immediately transferred to an m = ±1 (m = ^1) vortex in the signal (idler), leaving the pump 

vortex-less. 

^^ Later, in Sec. 



6.1 



in connection to the stability of multiply quantise vortices, we also describe 
vortices in the TOPO regime, where we follow the vortex dynamics not of the OPO like here, but 
of the extra population only. 
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signal, it has been possible to observe the vortex core slowly drifting, changing in 
shape and moving with different velocities. Note that, because these are metastable 
solutions, a minimum probe power is required for the polaritons to acquire enough 
angular momentum to be able to transfer it to the steady state. However, once the 
transfer is achieved, the probe power does not change significantly the duration and 
depth of the vortex in the steady state. 



5.2 Onset and dynamics of vortex-antivortex pairs 

There is an aspect that we have been neglecting in the discussion of the previous 
section on the occurrence of metastable vortex solutions in OPO triggered by an 
external LG probe. If the extension of the probe carrying a vortex with charge m ~ 
+ 1 is smaller than the size of the vortex-free OPO signal, continuity of the polariton 
wavefunction requires that necessarily an antivortex with charge m~—\ has to form 
at the edge of the probe (see Fig. \\A\ . Indeed, 'unintended' antivortices have been 
shown to appear in the signal at the edge of the imprinting vortex probe and we 
have explained in Ref. 1911 , both theoretically and via experiments, the origin of the 
deterministic behaviour of the antivortex onset and dynamics, i.e. where antivortices 
are more likely to appear in terms of the currents of the imprinting probe and the 
ones of the underlying OPO. 



5.2.1 Random phase between pump and probe 



As mentioned earlier in Sec. 5.1 single shot measurements would give a too low 
signal to noise ratio, therefore an average is performed over many pulsed experi- 
ments taken always for the same OPO conditions. What differs at each probe arrival 
is the random relative phase ^rdm between pump and probe, 

F (r , f ) = F„ (r, f ) + F„b (r , f )e'*^'"" , (26) 

with ^rdm uniformly distributed between and 2n. We simulate the dynamics 
of the vortex-free signal OPO (same conditions of Fig. 6 at fp — 1.2/?) follow- 
ing the arrival of a vortex probe ^A\ for 1000 realisations of ^rdm and then av- 
erage the complex wavefunctions over such realisations at fixed time and space. 

The steady state currents of the OPO signal before the arrival of the probe have a 
dominant component pointing leftwards and an equilibrium position where all cur- 
rents point inwards (bottom left part of the panel (a) in Fig.[T4|. In single shot sim- 



ulations of Fig. 14 d,f) (one realisation of the phase 'I'reim), we find that if the probe 
is positioned well inside the OPO signal, then the imprinting of a vortex m ~ +1 
(antivortex m = — 1) forces the system to generate, at the same time, an antivortex 
m = —\ (vortex m = + 1) at the edge of the probe. This is a consequence of the con- 
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Fig. 14 Profile and cuirents of the steady state OPO signal before the arrival of the probe (a) and 
associated interference fringes (d) — parameters for OPO are exactly the same as the ones of the 
inset 2 in Fig. 6 (/p = 1.2/*). Location of antivortices (dots (b)) and vortices (stars (c)) at the 
arrival of a vortex (stars (b)) or an antivortex (dots (c)) probe, for 1000 realisations of the random 
relative phase between pump and probe, <P,d,„- The size of dots in (b) (stars in (c)) is proportional to 
the number of times the antivortices (vortices) appear in that location. Panel (e) ((f)) shows single 
shot interference fringes relative to the plot in (b) ((c)). Contour-level lines in (b) and (c) represent 
the photonic disorder Vc{r)- The white circle represents the edge of the probe. From Ref. [91 J ask 
for copyright permission!. 



tinuity of the polariton wavefunctions: If the signal OPO phase is homogeneous and 
vortex-free before the arrival of the probe, then imposing a topological defect, i.e., 
a branch cut, on the signal phase at the probe core, requires the branch cut to termi- 
nate where the phase is not imposed by the probe any longer and has to continuously 
connect to the freely chosen OPO signal phase, i.e. at the edge of the probe. As re- 
peatedly mentioned in this review, OPO parametric scattering processes constrain 
the sum of signal and idler phases to the phase of the laser pump by 2(j)p = 0, + 0,. 
Thus, at the same positions where the V-AV pair appears in the signal, an AV-V pair 
appears in the idler, so that locally the phase constraint described above is satisfied. 
This agrees with the experiments in ID, though there only a single V (AV) in the 
signal (idler) could be detected, because the signal size was comparable to the probe 
one. 

Different relative phases <Prdm cause the antivortex (vortex) to appear in differ- 
ent locations around the vortex (antivortex) probe. However, on 1000 realisations of 
the random phase uniformly distributed between and 2n, we observe that the an- 
tivortices (vortices) are more likely to appear on positions where the current of the 
steady state OPO signal before the probe arrival and the probe current are opposite. 
For example, for the ;« = +1 (m = ~1) probe of Fig. 14 c) (Fig. 14 e)), the current 
constantly winds anti-clockwise (clockwise), therefore, comparing with the signal 



current of Fig. 14 a), the two are anti-parallel in the bottom right (top left) region 
on the probe edge, region where is very likely that an antivortex (vortex) is formed. 
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Note also that the onset of antivortices (vortices) privileges regions where the steady 
OPO signal has a minimal intensity. This agrees remarkably well with it has been 
recently measured experimentally, in Ref. ||9T| . 



5.2.2 Multi-shot averaged dynamics 

Crucially, via numerical simulations, we elucidate the reason why an experimental 
average over many shots allows detecting a vortex by direct visualisation in density 
and phase profiles. Recently, it has been suggested by stochastic simulations II2TI 
that vortices in non-resonantly pumped polariton condensates undergo a random 
motion which will hinder their direct detection, unless they are close to be pinned 
by the stationary disorder potential and thus follow a deterministic trajectory |7 1. In 
the case considered here of a superfluid generated by OPO, we can instead explain a 
deterministic dynamics of the V-AV pair in terms of the OPO steady state currents, 
which determine a unique trajectory for the pair, allowing their observation in multi- 
shot measurements. 



By averaging the 1000 images obtained at the probe arrival, e.g., in Fig. 14 c), 
neither the imprinted vortex nor the antivortex can be detected: Both phase singular- 
ities are washed away by averaging the differently positioned branch-cuts. However, 
the steady state signal currents push the V and AV, initially positioned in different 
locations, towards the same equilibrium position where all currents point inwards. 
Thus, exactly at the time where the probe is shined, on average there is no V-AV 
pair, after ^ 10 ps, both V and AV appear and last ^ 75 ps (see Ref. [91 1), till they 
eventually annihilate. 

It is interesting to note that it has been experimentally shown |91 1 that the onset of 
vortices in polariton superfluids does not require a LG imprinting beam, but instead 
vortex-antivortex pairs can be also generated when counter-propagating currents 
are imposed, similarly to what happens in normal (classical) fluids. In Ref. |91| a 
Gaussian pulsed beam has been shined either at rest with respect to the OPO signal, 
kph = kj ~ 0, or moving kp/, 7^ kj. While no vortex-antivortex pair appears in the 
first case, in the second, a vortex-antivortex pair appears on opposite sides of the 
probe edge. 



6 Stability of multiply quantised vortices 

The energy of a vortex is proportional to its quantum of circulation squared ll76ll . 
np-. Thus, ignoring interactions, a doubly charged m = 2 vortex, has higher energy 
than two single m = 1 vortices. However, including interactions between vortices, 
the energy of an m — 2 vortex turns to be the same as the energy of two m = 1 
interacting vortices close together The behaviour of doubly quantised vortices has 
been the subject of intensive research in the context of ultra-cold atomic gases. In 
particular, it has been established that the nature of the splitting is the dynamical 
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instability. Nevertheless, m — 2 vortices have been predicted to be stable for specific 
ranges of density and interaction strength 1921 |93l . though, so far, they have not been 
observed experimentally 1*941. As for single vortices, multiply quantised vortices can 
be however stabilised in multiply connected geometry. Indeed, stable pinned m = 2 
persistent vortices have been recently observed 1891 by using a toroidal pinning 
potential generated by an external optical plug, and demonstrated to split soon after 
the plug was removed. In this case, the presence of a plug beam at the vortex center 
can pin both m = 1 and m — 2 vortex states and stabilise them against respectively 
spiralling out of the condensate for m — 1 and splitting for m = 2. In other words, the 
external trap mechanically prevents the persistent flow to undergo any movement. 

In contrast to equilibrium superfluids, such as atomic gases, both stable and un- 
stable m = 2 vortices has been experimentally realised in polariton OPO superflu- 
ids Is). In this section we provide a theoretical explaination of the stability and 
splitting of doubly charged polariton vortices. As done previously, vortices in OPO 
are generated by an external pulsed probe ( p5| ). As such, we classify the response of 
the system to an m — 2 LG probe, depending whether the probe generates a TOPO 
state (and the vortex is only carried by the extra population but is not transferred 



into the OPO signal), as described in Sec. 6.1 or instead is transferred in the OPO 



signal (Sec. 6.2 1 



6.1 TOPO regime 

We first consider the TOPO regime (see Sec. HI), i.e., when the vortex propagates in- 
side the triggered probe and conjugate wave-packets. It has been found |8, 23] that, 
in the TOPO regime, m — 2 vortices are stable within their lifetime when triggered 
at small momenta kph (see Fig.flSlpanels (a) and (b)), while they split into two m = l 
vortices for large values of kph (see Fig.fTSlpanels (c) and (d)). This conclusion was 
reached both by experimental observations |8| and theoretical analysis |8 23 1. The 
numerical analysis shows that the crossover from non-splitting to splitting happens 
for the probe momenta where the LP dispersion deviates from the quadratic one (see 
Fig.|l2|. The two different cases are shown in Fig. 15 For kph = 0.2/im^', at short 



times, the probe propagates without changing shape and with little change in inten- 



sity (not shown in the Fig. 15 1, consistent with the linear dispersion of spectrum 
characterising this regime. However, at longer times the density of the triggered 
probe and conjugate states drops more then two orders of magnitude, the dispersion 
changes to the quadratic one and the wave-packet expands (panel (b)). A uniform 
expansion of the wave-packet leads to the decrease of the probe and the conjugate 
polariton densities and thus to an increase of the vortex core, but it does not cause 
the vortex to split. In contrast, for kph = 1 .4/xm^ ' , where the LP dispersion is not 
quadratic, the m = 2 vortex state splits into two m — 1 vortices shortly after the ar- 
rival of the probe (panel (d)). This behaviour can be understood by analysing the 
evolution of the system's excitation spectrum in time: the dispersion of the time- 
dependent TOPO evolves from LP (before the probe arrival) to linear (at early times 
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after the probe arrival, when the stimulated scattering is strong), and back to the 
LP at later times. For large kpi, the LP dispersion deviates strongly from quadratic 
(see Fig. 12 right panel). Wave-packets propagating with non-quadratic dispersion 



do not keep their shapes (as discussed in Sec. 4.1 1, and the simulations show that the 
distortion can be very pronounced in particular at later times of the evolution. The 
distortion during the early times of the propagation leads to the mechanical splitting 
of an OT = 2 vortex, analogous to the structural instability discussed in Ref. ||95]| . 
Additionally, as discussed in Ref. |8|, for small kph, within the quadratic part of the 
dispersion, the group velocity of the wave-packet carrying the vortex equals the ve- 
locity of the net super-current (given by kpt,) associated with phase variations. This 
is not however the case for larger kph, beyond the quadratic part of the dispersion. In 
this case, the propagating vortex feels a net current in its moving reference frame, 
which may provide additional mechanism for splitting. 




Fig. 15 Intensity of the TOPO probe profile \ilf^''{r,t)\ after tlie arrival (at f = 0) of an ;n = 2 
vortex pulsed probe l |25) with Opt, — 87/jm. Small kpi, is shown in panels (a) and (b), while a large 
kph in panels (c) and (d). While in the case (a,b) the m = 2 vortex does not split within its lifetime, 
in (c,d) the vortex splits soon after the probe arrives. The intensity scale in (b) is 200 times smaller 
then in (a) - signal expands as the density drops two orders of magnitude (see text). Adapted from 

da. 



6.2 OPO regime 



In contrast to the TOPO regime described above, it has been shown both experi- 
mentally and theoretically [8, 2X1 that m = 2 vortices that do get imprinted into the 
steady-state OPO signal are never stable and spUts into two m = 1 vortices almost 



immediately, even before the probe reaches its maximum intensity (see Fig. 16 1. By 



analysing the system's dispersion in different regimes, as well as the dynamics of 
currents visible in the simulations, we have been able to identify several causes for 
the splitting: Before the arrival of the probe, the steady-state OPO dispersion is flat 
around the pump, signal, and idler However, the triggering probe favours the signal 
and conjugate to lock and propagate with the same velocity v^^ . This behaviour 
corresponds to a linear dispersion. Further, once the vortex gets imprinted into the 
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Stationary OPO signal and idler, the system's dispersion changes back to be flat. The 
evolution of the dispersion between flat, linear and again flat leads to a complicated 
dynamics of both signal and idler (the transient period described in |22|), causing 
the structural instability and splitting of the m = 2 vortex during the transient time. 
Another reason for the structural instability and splitting are the non-uniform cur- 
rents (see Fig.[T6| present in the OPO signal caused by the interplay between spatial 
inhomogeneity, pump and decay, which the OPO vortex experiences in its reference 
frame. 




Fig. 16 Filtered signal profile and currents above threshold for OPO (fp = 1.12/,, ) for a top-hat 
pump with FWHM cjp = 35 /im at f = — 3 ps (first panel) before the arrival of an m = 2 probe. The 
doubly quantised vortex gets transferred from the probe into the OPO signal and splits into two 
m = 1 vortices even before the probe reaches it's maximum intensity at f = ps (second panel). In 
this simulation the vortices coexist for sometime (roughly 15 ps), then one gets expelled from the 
signal (fourth panel). Adapted from (23l . 



7 Vortices in other polariton fluids 

We do not pretend to give an exhaustive review of the broad field which studies 
vortices in polariton fluids in this chapter, where we have mostly focussed on the 
occurrence of vortices in polariton OPO superfluids. However, we would like at 
least to briefly mention what happens for polariton fluids other than OPO. 



Spontaneous vortices in trapped incoherently pumped polaritons 



For incoherently pumped polaritons, the presence of a harmonic trapping potential, 
can make the non-rotating solution unstable to the spontaneous formation of a vor- 
tex lattice 1 80 1 (this work has been generalised to include the effects of polarisation 
in Ref. |96|). In experiments, spontaneous vortices in incoherently pumped systems, 
have been observed in Refs. If5l l97ll and their existence explained in terms of pin- 
ning by the disorder present in the CdTe sample. Adding the polarisation degrees 
of freedom, can give rise to the appearance of half- vortices II98I I6I. Polariton vor- 
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tices have been also observed in cavity mesas lISTl . cross refer to Benoit 
Deveaud' s chapter? 

Vortex-antivortex pairs have been observed in the non-resonantly pumped ex- 
periments of Ref. 1 10], where the mechanism of V-AV generation is explained in 
terms of density fluctuations originating from the cw multi-mode pumping laser, 
while for a single-mode laser no V-AV pairs have been observed. In this sense, their 
motivation and interpretation is in terms of the BKT transition. The pair dynamics 
in the condensate has been studied in Ref. Il99l . cross refer to Yoshisha 
Yamamoto' s chapter? 

Finally, as mentioned previously in SecslSland 5.2.2 generation and detection of 



metastable vortices have been also recently discussed for polariton condensates gen- 
erated by toroidal non-resonantly pumping in Ref |21 J, where vortices have been 
seeded with an external LG probe. Interestingly, very recently in Ref. 1 100|, it has 
been observed an all-optical spontaneous pattern formation in a polariton conden- 
sate non-resonantly pumped with a ring geometry. 



Resonantly pumped-only polaritons 

In Ref. lfTn[T2 l. the hydrodynamic nucleation of V-AV pairs is studied by making 
collide the polariton fluid with a large defect. In particular, polaritons are resonantly 
(coherently) injected with a pulsed laser beam, creating a population in the pump 
state only. The focus and interest of these studies are the possibility of exploring 
quantum turbulence, the appearance of dissipation and drag above a critical velocity 
because of the nucleation of vortices in the wake of the obstacle. 
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